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Abstract 

We consider a two matrix model with gaussian interaction involving the term 
trABAB, which is quartic in angular variables. It describes a vertex model (in 



particular case - of F-model type) on the lattice of fluctuating geometry and is 
the simplest representative of the class of matrix models describing coupling to 



two-dimensional gravity of general vertex models. This class includes most of phys- 
ically interesting matrix models, such as lattice gauge theories and matrix models 
describing extrinsic curvature strings. We show that the system of loop (Schwinger- 
Dyson) equations of the model decouples in the planar limit and allows one to find 
closed equations for arbitrary correlators, including the ones involving angular vari- 
ables. This provides a solution of the model in the planar limit. We write down 
the equations for the two-point function and the eigenvalue density and sketch the 
calculation of perturbative corrections to the free case. 
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1 Introduction 

The matrix model approach to 2D gravity (see |I|, [| for a review) has proven to be 
very useful and illuminating, even providing a method of addressing the non-perturbative 
issues of string theory || . It suffers, however, a severe limitation, since all matrix models 
solved up to the present time correspond to strings embedded in spacetimes of unphysical 
dimensions c < 1. The models with c > 1 are investigated either perturbatively Q, 
or numerically 0. The same c = 1 - barrier appears in the continuum approach to 
non-critical strings too, making them inappropriate for description of systems like gauge 
theories |], [/J. Its physical origin is, at the present time, not very clear. There are 
indications that it reflects an intrinsic instability of random surfaces with respect to 
crumpling , is possibly connected to the existense of the tachionic mode of the bosonic 
string , and can be understood in terms of condensation of spikes ]7|, BL |9], [TTJ. These 



observations suggest that the corresponding phase transition may be controllable in a 



model of surfaces with the action depending on the extrinsic curvature |L_T], [TIJ. Being 
initially simpler, the matrix model formulation of such models could shed a new light on 
the problem and even suggest a way out of it, were it not for the technical obstacle to 
solving interesting multi-matrix models which consists in inapplicability of the large-N 
saddle point technique to the systems with the number of degrees of freedom of order N 2 
(N being the size of the matrix). The latter is not the case for one- matrix models, where 
the action does not depend on the angular part of the matrix and one deals with a system 
of N eigenvalues. For the next case, the simplest two-matrix model, Itzykson and Zuber 



were able to circumvent the problem by integrating over the angles explicitly [12]. The 
beautiful formula of Itzykson and Zuber allows for the solution of multi-matrix models 
with the tree-like structure of the target graph. A direct further generalisation of this 
approach to target spaces containing cycles would invoke calculation of the correlators 
of Itzykson-Zuber integral, for which only rather complicated expressions are known [13] 



(see, however, the interesting work of ref. [23). Given the importance of the matrix 
models, it is therefore of some interest to look for models which are tractable by methods 
not involving the elimination of angular variables, such as the method of loop (Schwinger- 
Dyson) equations (see [Lj| for a review). The class of such models includes, of course, the 



ones solvable by previous techniques, in which case solving the loop equations sometimes 
allows the calculation of previously unknown quantities, such as the correlators involving 



angular variables JT51 [HI ITS1 . The interesting question is whether it is possible to address 



in this way any new models. In the present letter we describe a step in this direction. 



We define a gaussian 2 matrix model describing coupling to gravity of a vertex model (in 
particular case - of F-model type, [|19[]). The action of this model contains a term trABAB 
which prevents one from using the approach of Itzykson and Zuber. However, the infinite 
system of planar loop equations of the model can be reduced to a closed set, thus giving 
an analytical, though not very explicit, solution of the model and allowing one to find, in 
principle, all planar correlators, including the ones involving the angular variables. This 
provides a new example of a situation where loop equations are more effective than the 
"standard" approach. 

The model we consider belongs to a very general class of matrix models describing 
vertex models coupled to gravity. This class includes most of the physically interesting 
matrix models, among them lattice gauge theories and the matrix models describing 
the extrinsic curvature strings. In view of the above, the latter property seems to be 
particularly interesting and is our primary motivation for introducing and considering 
this class of models. We comment on this connection in the next section, after having 
introduced our model which is, unfortunately, only the simplest representative of the class 
and corresponds to c = 1. We hope, however, that it shares some formal features with 
its more interesting relatives. It is interesting by itself, too, as a novel c = 1 two-matrix 
model. 

In the last section we use the orthogonal decomposition of the matrices (used recently 
in the works fl5| , [16], [17]]) to derive a system of loop equations written in the variables, 
allowing for its decoupling. We obtain a Fredholm equation of the second kind for the 
eigenvalue density and a more general non-singular nonlinear integral equation of Ham- 
merstein's type for the general two-point function < trA n B m >. This equations have 
unique solutions in terms of Fredholm series. Analogous equations exist for arbitrary pla- 
nar correlators of the model. The equations we obtain are rather complicated, but well 
suited for perturbation expansion. We sketch the calculation in the first order. 

2 A model 

Consider a system of two hermitian matrices with the partition function defined by 

Z= [dAdBe~ NS , (1) 



and the quadratic action 

S = tr-(A 2 + B 2 ) - IctrAB + -trABAB + -trA 2 B 2 . (2) 

2 y ' 2 2 y ' 

This action leads to a perturbation theory involving two types of propagators, 

P(A,A) = P{B 1 B) = -^— I , (3) 

P(AB) = -j-H-j, (4) 

and two types of vertices, with the couplings A and «;. One can imagine labeling, for 
each vertex of the given diagram, every A-type (B-type) leg by an outgoing (incoming) 
arrow. One then obtains a configuration of a vertex model, which is very similar to the 



familiar F version of the 6- vertex model [19[, which has a critical regime with c = 1. The 



difference in the present case is that the propagators of the type (|3]) allow for configurations 
with directions of the arrows not preserved along the link. This may be viewed as some 
deformation of the F-model. One can hope to recover the pure F-model continuing the 
final expressions to the point r = 0. Note, however, that this limit is singular from 
the point of view of the original integral (|l|). Alternatively, one can put from the very 
beginning c = 0. This will leave the diagramms, containing only the "wrong" propagators 
(0). Than one can change locally the direction of all arrows for half of the vertexes to 
show that there exists a configuration of the F-model corresponding to each diagram. For 
some diagrams this procedure will meet, however, global obstructions. Note also, that 
from the vertex model point of view the terms A 4 and B A would correspond, if added to 
the action, to deforming the 6- vertex model to a non-critical 8- vertex model |T^. A way 
to couple the F-model to gravity is to use the complex matrices which give rise naturally 
to the oriented propagators (see ||). This corresponds to taking c = and adding the 
quartic terms with the coefficients chosen to ensure the charge conservation. Such a model 
arises also as a reduction of the generalised Weingarten model and is solved in general, 
and, in the planar case, in the presence of a vortex term, in the works ^ZQ, ^TJ. [] 

The model (0) belongs to a very general class of matrix models described by the actions 
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S = -trl af3 <& a & + V^^W + ... (5) 



where V^V (with V^ = I) are invariant with respect to cyclic transposition of indices. 
In diagrammatic expansion of this model l a/3 ( = inverse of l a p) is the propagator of the 
matrix field $ and V^ with n > 2 represent the vertexes. From the point of view of 
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the lattice dual to the given diagram one has a statistical model with "spins" a living 
on the links of a triangulation of some surface, and the statistical weight obtained by 
prescribing the factor V^ to each face (triangle for n = 3) and l al3 to each pair of the 
links identified in forming the surface. This action was first written down by Bachas and 



Petropoulos in connection to topological models on a two-dimansional lattice p2| . The 
spins a take values in some set G. The familiar matrix models are the particular case with 
K»/3 7 ... ~ 'W^cry ■ ■ ■ ,(n > 2). One can take 

[$«]t = $« f (6) 



with the dagger on the left hand side denoting the hermitian conjugation and on the r.h.s. 
- some involution on G. For a^ = a one has a hermitian matrix model. 

Let us show that the matrix models of this type describe, for special choice of G, 
the extrinsic curvature strings. To see the idea take a to be a pair (x,n) where x is a 
coordinate in 3-dimensional space and n is a unit 3-vector, and arrange the couplings to 
have 

S= n J2 tr&(x,n 1 )&(x,n 2 )ln 1 n ll + '52trV(&(x,n)) 

2,ni,n2 x,n 

+ Y, ir$(a; 1 ,n)<l>(a: 2 ,n)<l>(x3,n) ( 5(n~x 12 xx 2 3)e- Ar ' ea(123) . (7) 

In the perturbation expansion of this model each link of the triangulated surface dual to 
a diagram is labeled by a position x in the embedding space of the string and there is a 
unit vector n ascribed to each face of the triangulation. The 6 -function in the cubic term 
ensures that n is a unit normal to the triangle (123), and one may take / ni ' n2 = exp(—nf 2 ) 
to have this mutual orientation-dependent factor associated with each neighbouring pair 
of triangles. The Area(123) denotes the area of the triangle (123), while the potential 
term represents the couplings which add to the intrinsic area of the surface leaving the 
extrinsic one unchanged, and thus describe the fluctuating intrinsic metrics (absent in the 



usual continuum space formulation of the extrinsic curvature strings [ 11 1 ) . 

This action, though illustrating how one can hope to obtain the matrix model de- 
scribition of extrinsic curvature string, is not the "final" correct version of the extrinsic 
curvature matrix model. One still has to modify it by taking the matrices to be complex 
and depending on some additional unit vector. For the discrete embedding space this 
corresponds to matrices living on the oriented links of the lattice and interacting along 
the plaquettes in the usual Wilson way; the index n will describes in this case the orien- 
tation of the plaquette. Identifying the matrices with different n one gets the Weingarten 



model p3) (see |21|). Generalization to higher dimensions in obvious. This connection 
of such matrix models to extrinsic curvature strings is our primary motivation for their 
introduction and investigation. 

Coming back to the model @, let us note that due to the presence of the interaction 
term trABAB the target graph of the model contains cycles preventing one from using the 
approach of Itzykson and Zuber. In spite of this, it is possible to proceed in calculating 
arbitrary correlators of the model. We will show below that the system of the planar loop 
equations of the model reduces to a finite set. 

3 Loop equations 

It turns out to be convenient to write the loop equations in terms of the variables defined 
via a decomposition 

M = Y^m aa P aa} (8) 

a 

where m aa are the eigenvalues and P aa the projectors on the corresponding eigenvector. 
The latter satisfy the conditions 

* act* a/3 OafiP aai ™ '*aa J-j / t *aa *-. \\3) 

a 

We will denote below the eigenvalues of A (resp. B) by x a (resp. yp) and the corresponding 
projectors as A a and Bp. Another useful notation is 

Pa 1 p 1 a 2 p 2 -a n p n = tr A ai B p x A a2 B p 2 ■ ■ ■ A an Bp n , (10) 

and also 

*nma 2 /32---a k p k " 2-~i X aiVpi on/3 1 a2/32---a k p k - {**) 

ati/3i 

To avoid the possible confusion let us keep in mind that the numbers in the sequence of in- 
dices of the objects like Pi2a 2 /32-a„0n will always denote the power of the matrix in the cor- 
responding position in the trace in flTT|), while the Greek letters will be used for numbering 
the projectors. The same will apply to the continuum notation P(l, 2, x 3 , y 4 , • • • , x n , y n ) 
with continious variables Xi (yi) taking the places of the discrete indexes oti [fii]- The 
variables (§) were recently used by Alfaro to solve the 2-matrix model and to write down 



the loop equations for D-dimensional matrix models Hl6| , |PT |. 



The equations we need follow from identities of type 

d 



J —(A^BfaA^Bfa ■ ■ ■ A an Bp n ))e- NS = 0. (12) 



In particular, one finds for n = 2 

rx ai — Cy/3 2 + —X ai Vi3 2 \Pai0ia202 ~ T7 Z^ ( "ai/3ia2/32 + P-yit3 1 a 2 P2J 

11 K 

— 1T7 z2 Pa 1 t3 1 P-y(3. 2 {^a 1 a2 ~ &ia 2 ) + Ml faPai Pitfall + ^Pa 1 f3 1 a2f3 2 12 = 0. (13) 

Summing this expression over a 2 and using (|9j) one obtains the relation for n — 1: 

k 1 1 K 

[rx Q - q//? + -x a y|]P a/3 - T7 S (^ + P-yfi) + typPapu + ~P a pi2 = 0. (14) 

An analogous set of equations is obtained by taking in (|T2"D the derivative -^ instead 
of J^. This new equations can be taken into account by requiring P ai f3 ia2 (3 2 --a n p„ to be 
invariant under the cyclic transposition of indexes. 

In the limit of N — > oo we introduce the usual rescaled variables Xi = x a , % = j?, 
yj = yp,j = % and assume that the eigenvalue density p{x) = ^- (the same for both 
matrices) has a finite support [—a, a]. One obtains, for instance, for eq . ([H|) : 

r k 2 r a u(z)dz lT ^. r a u(z)dz^. , 

[rx - cy + -xy 2 + 4- ^-^-}P{x,y) +4- ^^P(z,y) 

Z J —a Z — X J —a Z — X 

+XyP(x, y, 1, 1) + ^P(x, y, 1, 2) = 0, (15) 

with the density of eigenvalues and the correlators normalized by the conditions following 
from (191) 



p(x)dx = 1, (16) 

p(x)P(x,y)dx = 1, (17) 

f a p(x)P(x,y,l,l)dx = yP{l,y), etc. (18) 

J — a 

The symbol -f- denotes the principal value of the integral. The equation ( |T3"D is rewritten 
analogously as 

r K 2 f a u(z)dz r a u[z)dz . 

[rxi - cy 2 + -xiy 2 + -f \P(x 1 ,y h x 2 ,y 2 ) +f P(z,y 1 ,x 2 ,y 2 ) 

L J —a Z — X\ J— a Z — X\ 

u(z)dz 5(xi — x 2 ) 5(z — x 2 ) 1 



+P( Xl , yi )-L -±^P(z,y 2 )[ 

J —a Z — Xi 



Z-Xx ' p{xi) p(z) 

+XyP(x h yi, x 2 , y 2 , 1, 1) + -P(x x ,yx, x 2 , y 2 , 1, 2) = 0. (19) 



Taken for arbitrary n, (|i2l) thus provides an infinite set of linear relations between 



the correlators of the system. Our aim is to show, that eqs.fllSD and (|19"D reduce, in fact 



to a closed system. To see this, note that multiplying both sides of eq. ([15|) by p(x) and 
integrating over x one obtains an algebraic equation for the correlator P(x, 1) which gives 

P{x,\) = ^ . _ . (20) 

y ' r + (A + k)x 2 v ' 

A similar thing occures when one integrates over x\ in eq.(plj). One obtains the equation 

[r + \y lV 2 + ^{yl + y 2 2 )]P(l, yi, x 2 , 2/2) = 
J( yi -y2) _ r- dzp(^ [p( )p( j + p( )p( 5] (21) 

j0(yij ./-a 2 — X 

which expresses the 4-point correlator P(l,y 1 , £2,2/2) in terms of 2-point functions. From 
here one finds easily 

cu 2 
P(x, y, 1, 1) = r + / A+ x 2 P ( X ' ^) 

and 

[" iTv^T&T J ' ^lPMP(z,y) + P(x,y)PMl (23) 

J-a r + XQy + |((/ + y 2 ) J -a z — x 

Substituting these into eq. (|15]) gives a non-linear singular integral equation involving two 
unknown functions, p{x) and P(x, y). The second equation is provided by the symmetric- 
ity condition P(x,y) = P(y,x). 

Let us now rewrite this system resolving the singular part of the equations. We 
introduce the functions 

*(,) = f ^£, (24) 

J —a Z — X 

f a u(z)dz , , , 

<Kx;y) = / -±J—P(z,y), 25 

J-a Z — X 



and 



K 2 



G(x; y) = rx — cy — —xy + H(x). (26) 

The variable x in these functions takes values on the complex plane cut along the interval 
[—a, a], while y is real and belongs to [—a, a}. One has, for x G [—a, a], 

Disc x (j)(x; y) = 4>(x + i0; y) — <f)(x — i0; y) = 2irip(x)P(x, y) (27) 



and 

Cont x <f)(x; y) = <f>{x + tO; y) + <j>{x - tO; y) = 2 -f ^^P(z, y) . (28) 

J — a Z — X 

Using these relations, in (|15j), (0) and (|23|) one can easily show that the both sides of 
the equation flT5| ) are given by discontinuities of some functions. One obtains 

K CI I 

Disc x G(x; y)<j)(x; y) = -Disc x (X + -) , — T^O^ 2/) 

2 r + (A + K)y^ 

r a XCv + -C 2 

+Di S c x <j>(x,y) dCp(C) -TT^fe- ^ foO, (29) 

which implies that 

K" CI I 

G(x; y)(f)(x; y) + (A + -) — r^0(^, 1/) 

2 r + (A + K)y z 

r a XCv + -C 2 

-<j>(x, y) j_ a d(p(0 r+XC y + £ 2 + y2) <i>(x, = PoL(x), (30) 

where Pol.(x) is some polynomial of x. Taking x — ■> oo and using the known behaviour of 
if(x) and (f>(x,y) in this limit, one finds Pol.(x) = r — |y 2 and, after some algebra, (|30|) 
reaches its final form 

0(x; y)^(x; y) + [x - °^ - ]0(g; y) + 1 = (31) 

r + (A + K)y z 

with 

*»»^5' I( »^|m' (32) 

This is a non-linear integral equation of Hammerstein's type which defines the function 
p(x)(f>(x,y). Another equation, expressing the symmetricity of P(x,y) can be written as 
$(x,y) = $(y, x) with 

Hx,y)= f pp^cP(xX). (33) 

J-a Q-y 

Taking here the limit of y — > oo one obtains the simpler condition 

H(x) = f d(p(()<f>(x,z) 

dzp{z) (34) 



d(p(Qcf>(x,0 
r+{X+K)z 2 ' J-a r+ACy+f (C 2 +2/ 2 ) 



a # r+fA+K^ 2 + -/"-a 



Equations ( |3T]) and (^) allow one to find the eigenvalue density p(x) and the two-point 
function P(x, y), and provide the decoupling of the system flT2|). The existence of solutions 



for (pTD and fl33|) follows from the theorems of the theory of integral equations, but we 



do not know, at the present moment, how to find them explicitly. The same theorems 



say a lot about the properties of the solutions. Leaving the complete investigation of the 
system (|3l| , |33| ) to the future, let us have a brief look on its perturbative properties. 

Consider first the free case. Taking A = k = one finds from the symmetricity 

condition (J34|) 

p[z)dz 



-HJX 



or 



H, 



-a rx — cz + Hq{xY 
rx + H (x)^ 



c 



o 



H (x). 



c I r 

This functional equation is similar to the one found in different contexts by Matitsyn 



(35) 
(36) 



and by Boulatov ||25|| . It can be solved either by noticing directly that it implies 

rH (x) 2 + (r 2 — c 2 )xH (x) = const. (37) 

and finding, for x —> oo, the constant to be c 2 — r 2 , or by taking c <^C 1 and performing 
successive approximations. The solution is given by 



H (x) 



r 2 — c 2 



2r 



4r 



•x* 



ry> — ^^_ / ' — 



(38) 



and leads, of course, to Wigner's semi-circle distribution 



Po{x) 



r 2 — c 2 

2irr 



4r 



J r 



t * ^-i ^^_ Z^" 



The 2-point function is, in this limit, given by 



rx — cy — H (x) 



(39) 



(40) 



Coming back to eq. ([31]) note that in the limit of x — » oo the non-linear term in it is 
much smaller that the other two. Applying the successive approximations to 

1 1 



(f)(x;y) 



x 



cy 



r+(A+K)y 



.;; 



cy 



-(f>(x;y)K<j)(x;y), 



(41) 



r+(\+K)y 2 



one obtains <fi(x,y) as an expansion in 1 cy — . The first term is given by 



r+(A + K);/ 2 



<f>W(x;y) 

while the succesive ones contain p(x): 

1 



x — 



cy 



(42) 



<f>U(x;y) 



r+(X+K)y 2 



<M 



x — 



cy 



r+{\+n)y' 2 



X — 



cy 



<p w (x;y)K<p w (x;y), etc. (43) 



r-\-(\+K)y 2 
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Comparing fl43|) to the expansion 



(„) = -I - P ^ P(2 ^ 



ry nn£ rr"J 



(44) 



one finds the correlators P(n,y) in terms of p. In particular, 

d(p(() 



Pfrv) 



2 
C X 



+ 



r + (A + k)x 2 J- a r + \(y + §(C 2 + y 2 ) 
Note an exact relation following from (f43l) : 



(45) 



(x;y;c) = 0(x 



cy 



r + (A + K)y- 



r,y,c= o). 



(46) 



Finally, multiply both sides of eq. (|T^) by p(y) and integrate over y e [—a, a]. One 

obtains 

f p(z) 



f 

J —a 



1 



'—a 2 X 

Using (H) to find 



[rx - cP(x, 1) + kxP(x, 2) + AP(x, 1, 1, 1)] . 



P(x, 1,1,1) 



c 2 x 3 



+ 



^Cp(C)C 



[r + (A + /?)x 2 ] 2 i-ar + ACy + f (C 2 + I/ 2 



(47) 



(4J 



as well as (|45| ) and fl2"0|), one arrives at the equation containing p(x) as the only unknown: 

where 

f{x) = ~ rX - [r + (A +V] 2 + £ rfCp(C) r + A C x + + f (C 2 + * 2 ) ' ^ 

The same equation could, of course, be derived directly by integrating in ([[]) over one of 
the matrices. This would lead to an effective action involving 'non-local' interactions of 



the type (trA n ) m . Such matrix models were considered in the works of ref.[p6fl 
Once again, one can resolve the singular part of ( f49|) to reduce it to 

1 



H(x) = r 



2~7n~ 



y/a 2 — x 2 f a dz 



2m 



a Z — X y 



a* — z' 



-z + 



c 2 z 



dz 



a Z — X y i 



or — Z* J -a 



d(p(C) 



[r + (A + k)z 2 } 2 

X( + KZ 

r+AC^+f(C 2 + ^ 2 )' 



(51) 



The first integral in this equation can be calculated easily by rewriting it as a contour 
integral around the cut and inflating the contour to touch the singularities at z = x, 



±iJr/(\ + k) and oo. 



11 



Taking the discontinuity of both sides of (|51"D one obtains a Fredholm equation of 
second kind for p(x). One can then perform the succesive approximations once again to 
obtain the solution in terms of Fredholm series. Alternatively, one can use (|5T| ) directly 
to calculate H(x) as power series in A and k. Consider, as a simple example, the case of 
c = k = 0. One finds from (51), in the lowest order 



H(x) 



■>- 



x 



yx 2 



(52) 



where the coefficient p2 = J^ a dzp(z)z 2 , as well as a 2 , have to be determined from com- 
parison of ([52]) to the expansion 



/ X 1 P2 Pi 

P( x ) = -Z--^~ZK- 



One finds, after some algebra, 



P2 



x 


X 3 X 5 


a 2 

4' 


4 
r 


l + (-) 2 4 
p r z 



(53) 



(54) 

It is not difficult to write down the corrections of first two orders for general c. The 
formulae are, however, not very illuminating. 



4 Comments 

This concludes our preliminary investigation of the model (§). Evidently, there is still 
much to be done, as the progress achieved in the present work is mostly of technical 
nature. We feel however that it is of some interest, as a demonstration of usefulness of the 
loop equations for calculations of the correlators involving angular variables for the vertex 
model on the lattice of fluctuating geometry which was not considered previously. Our 
interest in the model (§) was triggered by the observation that it contains an interaction 
of the same type as the matrix models describing extrinsic curvature strings. The latter 
models can lead to important physical insights. Further work is in progress. 
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